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The quasigeostrophic model describes large scale and relatively slow fluid 
motion in geophysical flows. We investigate the quasigeostrophic model un- 
der random forcing and random boundary conditions. We first transform the 
model into a partial difi'erential equation with random coefflcients. Then we 
show that, under suitable conditions on the random forcing, random bound- 
ary conditions, viscosity, Ekman constant and Coriolis parameter, all quasi- 
geostrophic motion approach a unique stationary state exponentially fast. This 
stationary state corresponds to a unique invariant Dirac measure. 

1. Introduction 

The quasigeostrophic (QG) model is a simpHfied geophysical fluid model at asymp- 
totically high rotation rate or at small Rossby number. It is derived as an ap- 
proximation of the rotating shallow water equations by a conventional asymptotic 
expansion for small Rossby number II 3| . The lowest order approximation gives the 
barotropic QG equation, which is also the conservation law for the zero-th order 
potential vorticity. Warn et al. and Vallis |23 emphasize that this asymptotic 
expansion is generally secular for all but the simplest flows and propose a modified 
asymptotic method, which involves expanding only the fast modes. The barotropic 
QG equation also emerges at the lowest order in this modified expansion. 
Moreover, it has recently been shown El H] that quasigeostrophy is a valid 
approximation of the rotating shallow water equations in the limit of zero Rossby 
number, i.e., for asymptotically high rotation rate. The three-dimensional baroclinic 
quasigeostrophic flow model can be derived similarly; see, for example, [IHllTIf!^!??! 

E|. 

We consider the barotropic quasigeostrophic flow model ^] 
(1) AV^t + Ji^p, AtA) + p-ipx = z^AV - rA^ + W2 
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on a rectangle D ~ (0, 1) x (0, 1) C , where y, t) is the stream function, (3 > 
the meridional gradient of the Coriolis parameter, v > the viscous dissipation 
constant, r > the Ekman dissipation constant, W2 the noise due to wind forcing, 
and J the Jacobian operator, which is defined by J{f,g) = fxdy — fydx- Equation 
Q can be rewritten in terms of the vorticity q = /S.ip as 

(2) qt + J(V', 9) + Pi^x = i^Ag -rq + W2, 
which is usually supplemented with boundary conditions 

(3) ■4>{x,y,t) = on 51?, 

(4) -^qix^y^t) = Wi ondD, 

on 

along with an appropriate initial condition 

(5) q{x,y,0) ^ qo{x,y) , 

where n is the unit outward normal vector on the boundary dD, Wi is a temporally 
two-sided Wiener process with values in the function space U. The boundary condi- 
tion (O means no normal flow can pass through the boundary. The other boundary 
condition Q says that q = Aip has zero mean but with fluctuations, and thus 
might be called a random slip boundary condition 14 . As discussed in Pedlosky's 
book ( 14, , page 34), the boundary condition Aip = may be an appropriate slip 
boundary condition for the large scale quasigeostrophic motion. Boundary condi- 
tions for the quasigeostrophic model are not quite well understood, since this model 
describes large scale flows while boundary conditions also involve small scale mo- 
tions. For this reason we believe that, under the random media or random wind 
forcing conditions, a random slip condition may be more appropriate than the usual 
slip boundary conditions jj^ for the deterministic quasigeostrophic model. We also 
note that the Neumann form for the boundary condition on q is for mathematical 
convenience. 

In this article we treat the quasigeostrophic flow model with both random forcing 
and random boundary condition as a mathematical random dynamical system 
Our aim is to show that there exists a random steady state under a particular choice 
of parameter values. This random steady state is a statistically stationary solution 
towards which any other solution trajectory tends as t —> 00. Our mathematical 
approach is to formulate the random quasigeostrophic flow model as a stochastic 
evolution equation with structural similarities to the Navier-Stokes equation and 
then to show that it generates a random dynamical system for which there exists 
an attracting random fixed point. 
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2. Preliminaries 

Following Arnold ^ we will model noise in an abstract random dynamical system 
on a state space H hy a metric dynamical system {9, V) on a probability space V 
= (il, J-", P). A metric dynamical system consists of a group {0t}tem of operators Ot 
: 1^ fl, i.e., satisfying 

6*0 idn, 9s+t ^OsoOt for all s, i e R, 

such that the mapping (i, uj) ^ OtUJ from R® into is (S(R) x J^)-measurable 
and the probability measure P is ergodic (hence invariant) with respect to the flow 6. 

A random dynamical system consists of a metric dynamical system (0, V) and a 
cocycle mapping (p: RxflxHi-^H, i.e., satisfying 

(p{Q, uj, •) = id^f, ipis + t, Lu, •) — ip{s, Ot,-) o (p(t, oj, ■) for all s,t € R"*", 

that is {B{M.) ®J^®B{H), S(i?))-measurable. The mapping describes the dynam- 
ics of the system in the state space H , which will be a separable Hilbert space with 
inner product (•, •) and norm | • | = Vvv) iii this article. 

Let A" be a random variable defined on (9, V) with values in H . By the invariance 
of P and the measurability of the mapping (i, oj) i— > X{6tijj) from R x f2 into iJ is a 
measurable stationary stochastic process. We will restrict attention here to random 
variables generating stationary processes that satisfy certain growth conditions. A 
i?-valued random variable X is said to be tempered with respect to a ^-invariant set 
O' (of full P-measure) if the mapping 1 |A(6'j[jj)| grows at most subexponentially 
as i — > ±cx), i.e., for which 

\og+\x{etu)\ 

t^±QO \t\ 

for UJ G il' . Note that the only alternative to this when X is not tempered is 

log+|A(g,a>)| 
limsup j— j = -l-cio. 

t-+±oo \t\ 

A random fixed point of a random dynamical system is a _ff -valued random 
variable X* for which 

(6) ip{t,uj,X*{uj)) = X*{etLu) foraU<eR+ 

and for all a; in a ^-invariant set of full P-measure. The dynamics thus follows a 
stationary regime if we start in X*{u!). In particular, the probability distribution 
of these states is independent of t. 

The following theorem is a special case of a random fixed point theoremdue toSchmal- 

fuss [ni. 
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Theorem 2.1. Let tp be a random dynamical system with a separable Banach space 
{H, I'D as its state space. Assume that the mapping x i—np{t,LU,x) is continuous for 
every t > and lu E fl. In addition, let X — {X[ijj)}i^(z^ be a closed random set in 
H such that 

(7) E sup i,g k(l.^.M-^(l,^,MI ^0 

hi^h-2ex{i.o) \hi-h2\ 

and let the real valued random variable defined by u ^ ^^Vhex(u]) 1^1 tempered. 
Moreover, assume that 



\ip{t,uj,hi) - ip{t,u,h2)\ 

(8) sup sup 



te[o,i] hi^h2ex(uj) \hi - 

is tempered with respect to {9n}nei.- 

Then there exists a 9-invariant set of full ¥~measure fl' and a random variable X* 
satisfying |^ on f2'. Futhermore, X* is exponentially attracting, i.e., 

\un\ip{t,uj,X{uj))-X*{etUj)\^Q, a.s. 

t — >oo 

exponentially fast for any measurable selection X of X. 

2.1. Linear stochastic evolution equations. In the following we will consider 
the motion relative to a spatially constant flow. By the particular structure of the 
coefficients this spatially constant flow can by calculated separately. To find this 
spatially constant flow we have to solve a simpler equation. Henceforth we take for 
H the space of square integrable functions L2{D) on the rectangle D — (0, 1) x (0, 1) 
in fulfilling f{x)dD = and denote its norm by | • |. We then define V with 
norm |j • || to be the Sobolev space of functions contained in H with generalized 
derivatives of first order belonging to L2{D), and define W2{D) to be the space 
of functions with first and second generalized derivatives belonging to L2{D). In 
addition, we denote by J7 = L2(dD) a boundary space associated with square inte- 
grable functions on the boundary of D. To use the results of DaPrato and Zabczyk 
[T5] Chapter 13 we use that this perturbation is only defined on one side of the 
rectangle, say {0} x (0,1). So U consists of L2-functions on dD which are zero 
outside of {0} x (0, 1) with zero average. However, generalizations are possible. 

Let A be the Laplacian operator on D. The boundary value problem 

d 

^v^u = f, -7T—U — g on dD 
on 

with / G iJ, g G L2{U) and v > Q, has a unique solution u ~ G{f,g). The solution 
operator G : H x U 14^1(13) is a bounded hnear operator, i.e., there exists a 
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constant cg such that 

\\G{f,g)\\wiiD)<CG{\f\ + \\9\\u). 

Similarly we can consider the same equation but with the homogeneous Dirichlet 
boundary condition 

(9) -Au = f, u\eD = 0. 

The solution operator G for this boundary value problem satisfies 

\\G{f)\\wiiD)<CG\f\. 

In the following wc will denote by Cc^the constant which estimates G{f)x with 

respect to _L2(-D)-norm of /. 

In order to introduce a white noise on the boundary dD, we consider a temporally 
two-sided Wiener process Wi with values in U and denote by {!Ft}tes. the filtration 
corresponding to this Wiener process, where, roughly speaking, J^t is generated by 
the increments of the noise sample paths between times — oo and t. The associated 
Wiener measure P is defined on the cr-algebra of the canonical sample space 51 that 
consists of continuous functions cj from K into a phase space of the noise satisfying 
w(0) = 0. We will assume that the covariance operator with respect to this 
measure satisfies tiuQ^^ < oo and define the Wiener shift by 

Wi{9tu;,-) = Wi{uj,t + -)-Wiiu;,t). 

The measure P is ergodic with respect to the flow 9 of the metric dynamical system 

formed by the Wiener shift. 

Since the solution operator G(0, ■) above is a linear and bounded operator, the 
process t i— > Wi{co,t) := G{0,Wi{u},t)) also defines a Wiener process with trajec- 
tories in the Sobolev space Wf(£>) for which the bounded covariance operator Q 
:= G o o G* has a finite trace with respect to Wf (-D). 



The operator —//A with the vanishing Nemnann boundary condition can be ex- 
tended to an operator A defined on D{A) — VF^ (-D) with the vanishing Neumann 
boundary condition. The space H has a complete orthonormal base, consisting of 
eigenvectors ei, 62, • • • with corresponding eigenvalues Ai < A2 < • • • for the oper- 
ator A. By the particular choice of H we know that Ai > and that A is coercive. 
Let the semigroup {S{t)}t>o on H be the solution operator (indexed by t) of the 
initial-boundary value problem 

du du 

— - uAu = 0, m(0) =^06-^, = on dD. 

dt on 

with u > Q. This semigroup has the generator —A. 
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We consider an expression of the form 



z{uj,t) = S{t)zo+ ASit~T)dWiiuj,T), 



(10) 



as the solution of the hnear stochastic evolution equation 



(11) |-A. = 0, 



see (jni Section 13.2). The expression HlOfl is meaningful if, for example, 

[ \\AS{t - t)GQ^' i\\l2m H)dT < oo, for aU < f < +cx). 
Jo 

Since by the invariance of the increments of the Wiener process, we then have 



E 



A 



S{~T)dWl 



E 



A I SiT)dWi 



< oo. 



The random variable zwi defined by 

zwA^)-^a( S{-T)dWl{uj,T) 



is thus well defined and has finite second moment with respect to the norm of H . 
Moreover, we have formally 



S{t)zwA^) + A I S{t-T)dWi{uj,T) 



S{t)A S{-T)dWl{uJ,T)+ A S{t^T)dWl{uJ,T) 



A 



A 



A 



— oo 

t 



Sit- T)dWl {lO,t)+A S{t- T)dWl (w, t) 

Jo 

S{t-T)dWiiuj,T) 

S{-r)dWi{etU,T) = zwA^tio), 



so the stationary process 

t >^ z{etuj) ^ A f S{-T)dWi{9tuj,T), teR 

J —oo 

solves the boundary value problem (|ll(l . 

Since Elz^/i |^ is finite, we can apply the Burkholder inequality to obtain 

E sup jz^Vj (6't(jj)p < oo, 
te[04] 

and it then follows from the Birkhoff Ergodic Theorem ^S] that 



lim 

i— >±oo 
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on a ^-invariant subset of fl of full P-measure. Hence 

lim hldM^.o 

t-+±oo t 

on a 0-invariant subset of Q of full P-measure, i.e., \zwi \ is tempered. Note that 
similar techniques can be used to show that is defined on a 0-invariant set of 
full measure. 

Finally, equations for the generalized spatial derivatives of zwi can be investigated 
if we suppose that the covariance is sufficiently regular. Conditions are for- 
mulated in DaPrato and Zabczyk Theorem 13.3.1. In particular, Vzwi is well 
defined and tempered. 

3. Transformation of the quasigeostrophic equation 

We return to the QG vorticity equation in which we now write A^p for the 
vorticity. That is, we consider 

(12) ^ + J(V^, AV-) + Pi'^ = i^A V - r At/; + W2 {to, t), 

at 

with a non zero boundary condition 
d 

(13) —A^{t) = Wi{uj,t), V = ondD, 
on 

that involves a white noise Wi on the boundary dD as described in the previous 
section. In addition, the wind forcing white noise W2 is based on a temporally two- 
sided noise adapted Wiener process W2 with values in V and covariance such 
that tryQ^^ < 00. In particular, Wi and W2 are assumed to be independent. 
We can now define a metric dynamical system with the properties of our white 
noise terms. For Q, we choose an appropriate subset of the function space Co(R, U) 
X Co(IR., V^) with the usual Borel cr-algebra of a Frechet space, i.e., an element 
w is a continuous path from R into {U^V) with = 0. Then we take P = 

pWi ^ pW2 to be the product measure of the Wiener measures corresponding to 
Wi and which is ergodic since both and are ergodic. The flow 6 on 
is defined in terms of shift operators 9 applied to the sample paths of Wi and W2 ■ 

The above QG equation has structural similarities to equations of Navier-Stokes 
type. To be able to adapt well known results of such equations, we need to replace 
these boundary conditions by zero boundary conditions, which is possible with par- 
ticular types of stationary transformations; see Crauel and Flandoli [S] or Brannan, 
Duan and Wanner or in a more general context Keller and Schmalfuss or 
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Imkeller and Schmalfuss ^]|. In particular, we transform l|12|) into 

(14) + J(G(u),u) + PG(u)x = iyAu~ru + W2(uj,t) 
at 

d 

—uit)\^Wi{uj,t) ondD 
on 

where G is the solution operator of the boundary value problem (of. 

A-^ = u, ^p\dD = 0, 

i.e., with the solution i/; — G{u). We consider equation (|14|l as an evolution equation 
on the triple V C H C V where V is the dual space of V. 

The properties of the nonlinear term of equation (|14|) follow from those of the 
bihnear operator B : L2{D) x W2{D) V defined by 

(15) B{vi,V2) = J{G{vx),V2). 
Lemma 3.1. B is a well defined, continuous operator and 

i) {B{V1,V2),V3) ^ -{B{vi,V3),V2), 

ii) {B{vi,V2),V2) = 0. 

for vi e L2iD), V2 e W^{D), V3&V 

Proof. There exist positive constants c, c' and cb such that for any vi e L2{D), 
V2 e W2iD) C L4{D) and V3 e V C Li{D) we have 

\{B{vi,V2),V3)\ = {G{vi)xV2yV3 - G{vi)yV2xV3)dD 

Jd 

< c\NGiv,)\\L,iD)Nv2\\\v3\\L,iD) 

< C'\\\/G{V,)\\W,HD)\^V2\\\V3\\ 

< c_B |ui I |Vw2 1 ll-us II, 

which implies that B is well defined and continuous. 
Property i) follows from the integration by parts formula: 



G{vi)xV2yV3dD - G{vi)yV2xV3dD 

D JD 

- / G{vi)xV3yV2dD + / G{vi)yV3xV2dD 
Jd Jd 



G{vi)xV2V3Cos{n,y)dS — / G{vi)yV2V3COs{n,x)dS 

dD JdD 
= -{B{vi,V3),V2) 

because the boundary integrals are zero. Indeed, for two sides of dD these integrals 
are zero by the orthogonality of n and the direction of the derivative. For the other 
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both sides the integrals are also zero. For example, for the first integral we have by 
the properties of G 



Remark 3.2. If vi e V , which one can assume to be the solution of 1)14(1 . then 
we can similarly get that B{vi, vi) G C L2{D). This shows that we can split up 
the solution of the original equation into a special constant part plus the remaining 
part. Similarly we get G{vi)x G H. 

Equation (|14|l is similar to the equations of the Navier-Stokes type. Indeed, the 
Laplace operator term in ((14|l is also present in the Navier-Stokes equations (see 
Temam j2II|)j while the bilinear operator B defined by p5|l has similar properties 
(actually, a bit stronger) to the bilinear operator defining the nonlinearity of the 2- 
dimensional Navier-Stokes equations. It thus follows from the general theory of the 
stochastic Navier-Stokes equation that H14|l has a unique solution, see for instance 
Schmalfuss jl6| . The linear terms ru and (3G{u)x appearing in H14|l but not in the 
Navier-Stokes equation are not essential for a proof of existence and uniqueness. 
See Brannan, Duan and Wanner |4j for another proof of existence and uniqueness 
based on mild solutions. 



We now transform the stochastic evolution equation into a random evolution 
equation va V <Z H ^ V' , i.e., with stationary random coefficients rather than 
white noise driving or boundary terms. This will make it easier to find a forward 
invariant random set on which we can verify an appropriate Lipschitz condition. 
We introduce the random variable 



which we note without proof is a tempered random variable on a ^-invariant set 
of full measure. We also assume that W2 (hence Q^^) is sufficiently regular such 
that the Neumann boundary condition is fulfilled. Since 214^2 fulfills the Neumann 
boundary condition there is no influence to the boundary condition of (|14|l . 
We consider the random evolution equation 



G(fi)[x,l] = G(«i)[x,0] = 0, 
hence G{vi)j,[x,l\ = G{vi)x[x,Q\ = . 

Property ii) is a consequence of the antisymmetric nature of property i). 



□ 



4. The stationary solution 




—z{t, uo) + B{z, z) + Az + l3 G{z)^ + rz 



dt 



(16) 



= - B{z, zwA^t^) + zw2(6'tw)) - B{zw^ {Otto) + z^/^ (6*4^), z) 
- B{zwi (Otuj) + zw2 (^'t^), zwi (^t^^) + (OfUj)) 
~ PG{zwi (Oti^) + ZW2 {Ott^))x ~ r{zwt (Otuj) + ZW2 {Ot^j)) 
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with z{0) = zo <E H. 



Lemma 4.1. The random evolution equation has a unique solution for any 
initial condition zq G H and this solution defines a random dynamical system with 
respect to the metric dynamical system 9 introduced in Section 3 for which the 
associated cocycle mapping is defined by (t,uj,zo) ^ z{t,uj,ZQ). 

For the proof of this lemma we can use the fact that equation is quite similar to 
the Navier-Stokes equation. Although some linear terms are also present, similar a 
priori estimates can be obtained to those in Temam f|2()j. Chapter III) or Benssou- 
san and Temam 2 . Because of the properties of the operator B introduced in the 
previous section. Moreover, the random terms appearing inside the coefficients of 
equation 1)16(1 are given by stationary processes, so we obtain a random dynamical 
system, see Arnold (P, page 58). 



Remark 4.2. To see that t i-^ z{t,uj, zq) is continuous for any zq G H and oj € 
we can use Lemma IIL1.2 in since the solution of equation H16|l satisfies 

\\z{T,uj,ZQ)\\'^dT < oo 



for any zq <G H. Indeed, by the chain rule, 

\z{t)\'' + 2u f \\z{r)fdT < \zo\^ + 2{(3cG,.-r) f \z{T)\^dT 
Jo Jo 

+ 2CB I \z{T)\\\z{T)\\\S/zwA^r^)+^ZW2{0r^)\dT 

to 

||/(r)||_i||z(r)||dr, 

where / consists of all the terms in ((16|l that do not contain z. Then, using 

2\z\\\z\\\Wzw,+^zw^dT < '^\\z\\^ + -\z\^\\/zw,+Vzw,\\ 

2 V 

2||/ll-iPII < ^PlP + -||/f-i, 
2 V 

the asserted estimate follows by an application of the Gronwall inequality. Moreover, 
by the properties of the operators A and i?, we also have 

||z(t, tj,zo)||y'C^'^ < Zq ^ H. 
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We now define the random isomorphism i{Lo) : H ^ H hy 
i{ijj)a = -{zwiii^) + ziv2(^)) + o.^ 
for which the inverse isomorphism is given by 

z~^(Lj)a — {zwi{t^) + zw2{^)) + a. 
Note that the random variable i{uj)a{ijj) is tempered for any tempered a{oj). 
Lemma 4.3. Let z{-,u,zq) be the solution of J^lh]) . Then the process 

u{t, u), Mo) = i^^{9tuj) o z{t, uj, i{Lij) o uo) 
solves il4^ . In particular, u satisfies the boundary conditions i) j^] . 
Proof. The assertion follows by replacing z by u — zw-i{Ot^) ~ zwii^f-^)- D 

We will now check in the following Lemmata that the assumptions of the random 
fixed point theorem 12. II are satisfied. First, we show that there exists a tempered 
random set X(lu) of (single valued) random variables that will be mapped into itself. 

Lemma 4.4. Let X(uj) be the ball _B(0,p(a;)) in H with center zero and J-q- 
measurable radius 



p{uj) 



J exp (^(Aii/ - 2/3Cg,x + 2r)T 
+ \WzwdOsUj) + '^zw,{dsi^))\^ds^ ■R{9rUj)dT 



whe 



V 



and suppose that 



3c2 



where Xi > is the first eigenvalue of the operator A. Then the random set X is 
forward invariant, i.e., 

z{t,uj,X{uj)) C X{0tLu), t>0. 
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Proof. We have to estimate jzp for which we need the foUowing relations that are 
a consequence of Lemma 13.11 

2{B{z, z), z) = 0, 2{Az, z) = 2i/||zf > iy\\z\\^ + iyXi\z\^, 
2l3{G{z)^,z) < 2PCg,M^. 2r(z, z) = 2r\z\' , 

2\{B{z,zwA0tUj) + zw2{St^^)).z)\ < 2cB\z\\VzwAStUj) + '^zwAStUj))\\\z\\ 
< ^iVzwAOtLu) + VzwM^))\'\zf + ^INf , 

2{B{zw^i9tUj) + zwA^tL^)), z), z) ^ 0, 

2\{B{zwi{0tUj) + zw2{(^tuj),zwi{0tuj) + zvf^ (6*4^)), z) | 

< —\zwAStio) + zw,iSti^)\^\'^zwAStio) +VzwAOtu;)\^ + -\\z\\\ 



2\{PG{zwiidti^) + zw2{dtuj))x - r{zwi{dt^) + zi^^ (6'fti;), z) | 
< 3 ^- |zw,(0fw) + zwAdti^)\ + -ir\ 



It can be shown by a comparison argument that |z(t, zo)P is bounded by a 
solution of the afhne random differential equation: 

^ + (All. - 2I3Cg,x + 2r- zw^iOtLo) + VzM/.(^^t^))P))C = R{dtu) 

(17) '^^ 

m = \z,\\ 

for which the solution is given by a variation of constant formula. A direct calcula- 
tion verifies that t — > p^iOtto) is a solution of equation p7|l with initial value C(0) 
= p^{uj), which means that is a random fixed point of (|17|l . It thus follows that 
z{t,LU,ZQ) g X(9tuj) whenever zq G X(lu). □ 

We note that the random variable p is tempered (see ^Tj , page 110), so any selector 
contained in X is also tempered. 

It remains to check that the contraction condition Q of the random fixed point 
theorem holds. 

Lemma 4.5. Suppose that 

- J/Ai + 2PCg,x -2r+ — VzH'i + Vz^^ P + 



^^^^ 2r2 1 r2 2F /? 



, (1 + 2A^ Cg x/3)Ep' + — Ep^ + -^E| Vzwi + Vz^-, 1^ + < 0, 

where R was defined in Lemma \4.4\ Then the contraction condition ^ is fulfilled. 
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Proof. It follows immediately from Lemma ll-i. II with zi, Z2 ^ V that 

\{B{zi,Zi) - B{Z2,Z2),Z1 - Z2)\ = |(B(zi,Zi),Z2> + (-B(z2,Zi),Z2)| 

< \{B{zi - Z2,Zi),Zi - Z2)\ 

< Cb\zi - Z2\\\zi\\\\zi - Z2\\. 

Set zi = z{t,LU,hi), Z2 = z{t,LU,h2) £ H and 6z = zi — Z2- By the chain rule we 
obtain 

1^(1)12 + 21//" \\6zfdT < \hi-h2\^ + f {2Cg^x0\Sz\'^ ~2r\Sz\^ 
Jo Jo 

+2cb\Vzw, {0s^) + Vzw, {esUj))\\\6z\\\Sz\ 
+2cB\\5z\\\\zi\\\5z\)ds 

for hi, /i2 S X{uj). It then follows from 

2cb\VzwASsUj) + VzwM^))\\\Sz\\\Sz\ < M|Vzh^^ + VzwJ'|<5zp + ^||<5zf 

2cb\\6z\\\\zi\\\Sz\ < '^WSzf + ^llzif > : l\z\^ 

2 V 



that 



1 z(l,tJ, /li) - 2;(l,LJ,/l2)p^ 

-E sup log —2 

^ /ii5^h2eA'(cj) ^ "2| 



is less than or equal to the expression on the left hand side of inequality l|18|l. Here 
we have used the fact that 



2r) / p{B,ujf 
Jo 



V 



2 /g R{6rUj)dl 



\\Izw^{Ori^) + ^zwAdrUj)\ dr- 



for zo G -^(w). □ 

Remark 4.6. Note that the assumption of Lemma [4. 51 is sufficient for H4.4|l . 

The crucial point for the assumptions of the last lemma is to show at least for large 
V and for small iiuQ^^ and try^^Q^ that the random variables and p** have finite 
and sufficiently small expectations. In addition, the finiteness of the expectation of 
these random variables ensures that |(SJ) is satisfied. 

Lemma 4.7. The expectations of and p"^ are sufficiently small when v is suffi- 
ciently large and IyuQ^^ and tryQ^^ are sufficiently small. 
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We give only a brief comment on the proof of this very technical lemma. The 
essential ingredient is that and Vzvy^ are Gaussian processes, so R has finite 
moments of arbitrary order and 

(19) Ee^p(^aiJ^ \Vzw^\^dT^ < e"^*, 

and similarly for Vzwj- The constant ai > depends of the data v, Xi, ■ ■ ■ of the 
problem. The assertion of the Lemma follows if a2 is sufficiently small, which can 
be controlled by the traces of and . Finally, to obtain a finite dimensional 
version of the estimate (I19II we refer to Hasminskii (|H], page 37, Lemma 7.2), where 
we need the main assumption 

for an appropriate constant c. The variable can be handled similarly. 
Summarising, we have 

Theorem 4.8. Suppose that the assumption of Lemma \4.5\ is satisfied and let z* 
be the random fixed point of the (transformed) random dynamical system generated 
by Then there exists a random fixed point for that attracts the states of 

the phase space exponentially fast. 

Indeed, the random variable u that generates an exponentially stable stationary 
solution is given by 

u*{uj) = z*{uj) + zw2{^) + zwA^)- 
5. Discussions 

We have shown that, under suitable conditions on the random forcing, random 
boundary conditions, viscosity, Ekman constant and Coriolis parameter, all quasi- 
geostrophic motion approach a unique stationary state exponentially fast as time 
goes to infinity. In deterministic systems a high level of stability is obtained when 
there is an exponential attractor which attracts trajectories exponentially fast. In 
some situations this attractor is a single point (point attractor) which describes 
the laminar behavior of the flow. We are looking for such stability in the case of 
quasigeostrophic fluid motion under random perturbations. In particular, we flnd 
a random attractor which is defined by a single random variable. This random 
variable attracts all other quasigeostrophic motion exponentially fast. This random 
variable corresponds to a unique invariant measure, which is the Dirac measure 
with the random variable as the random mass point; see [T]. The corresponding 
stationary Markov measure is the expectation of this random Dirac measure. 



STABILITY OF THE STOCHASTIC QUASIGEOSTROPHIC EQUATION 



15 



References 

[1] L. Arnold. Random Dynamical Systems. Springer-Verlag, Berlin, 1998. 

[2] A. Bensoussan and R. Temam. Equations stochastiques du type Navier-Stokes. J. Fund. 

Anal, 13:195-222, 1973. 
[3] A. J. Bourgeois and J. T. Beale. Validity of the quasigeostrophic model for large-scale flow 

in the atmosphere and ocean. SIAM J. Math. Anal, 25:1023-1068, 1994. 
[4] J. R. Brannan, J. Duan, and T. Wanner. Dissipative quasigeostrophic dynamics under random 

forcing. J. Math. Anal. Appl., 228:221-233, 1998. 
[5] H. Crauel and F. Flandoli. Attractors for random dynamical systems. Prob. Th. Rel. Fields, 

100: 365-393, 1994. 

[6] B. Desjardins and E. Grenier. Derivation of quasigeostrophic potential vorticity equations, to 

appear in Adv. Diff. Eqns., 1998. 
[7] P. F. Embid and A. J. Majda. Averaging over fast gravity waves for geophysical flows with 

arbitrary potential vorticity. Comm. PDEs, 21:619-658, 1996. 
[8] R. Z. Hasrninskii. Stochastic Stability of Differential Equations. Sijthoff & Nordhoff, Alplien 
aan den Rijn, 1980. 

[9] D. D. Holm. Hamiltonian formulation of the bajroclinic quasigeostrophic fluid equations, Phys. 
Fluids 29 (1986), 7-8. 

[10] P. Imkeller and B. Schmalfuss. The conjugacy of stochastic and random differential equations 

and the existence of global attractors. Submitted, 1999. 
[11] H. Keller and B. Schmalfuss. Attractors for stochastic hyperbolic equations via transformation 
into random equations. Institut fiir Dynamische Systeme, Universitat Bremen, Report 448, 
1999. 

[12] P. Miiller. Stochastic forcing of quasi-geostrophic eddies. In P. Miiller R. J. Adler and B. Ro- 
zovskii, editors, Stochastic Modelling in Physical Oceanography, pages 381-396. Birkhauser, 
Basel, 1996. 

[13] J. Pcdlosky. Geophysical Fluid Dynamics. Springer-Verlag, Berlin, 1987. 
[14] J. Pedlosky. Ocean Circulation Theory. Springer-Verlag, Berlin, 1996. 

[15] G. Da Prato and J. Zabczyk. Ergodicity for Infinite Dimensional Systems. Gambridge Uni- 
versity Press, Cambridge, 1996. 

[16] B. Schmalfuss. Qualitative properties of the stochastic Navier Stokes equation. Nonlinear 
Analysis TMA, 28:1545-1563, 1997. 

[17] B. Schmalfuss. A random fixed point theorem and the random graph transformation. J. Math. 
Anal. Applns., 225(1):91-113, 1998. 

[18] S. H. Schochet. Singular limits in bounded domains for quasilinear symmetric hyperbolic 
systems having a vorticity equation. J. Diff. Eqns., 68:400-428, 1987. 

[19] T. G. Shepherd, T. Warn, O. Bokhove and G. K. Vallis. Rossby number expansions, slaving 
principle and balance dynamics. Quart. J. Roy. Met. Soc, 121:723-739, 1995. 

[20] R. Temam. Navier-Stokes Equation-Theory and Numerical Analysis. North-Holland, Ams- 
terdam, 1979. 

[21] G. K. Vallis. Potential vorticity inversion and balanced equations of motion for rotating and 
stratified flows. QuaH. J. Roy. Met. Soc, 122:291-322, 1996. 



16 



JINQIAO DUAN, PETER E. KLOEDEN, AND BJORN SCHMALFUSS 



(Jinqiao Duan) Department of Applied Mathematics, Illinois Institute of Technology, 
Chicago, IL 60616, USA 

E-mail address, Jinqiao Duan: duanSmath. clemson.edu 

(Peter E. Kloeden) Department of Mathematics, Johann Wolfgang Goethe University, 
D-60054 Frankfurt am Main, Germany 

E-mail address, P. E. Kloeden: kloeden@math.uni-frankfurt.de 

(Bjorn Schmalfuss) Department of Applied Sciences, University of Technology and Ap- 
plied Sciences, Geusaer Strasse, D-06217 Merseburg, Germany, 
E-mail address, Bjorn Schmalfuss: schmalfuss9in.fh-merseburg.de 



